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The mean field (MF) approximation for the pion matter, being equivalent to the Leading ChPT 
order, involves no dynamical loops and, if self-consistent, produces finite renormalizations only. 
The weight factor of the Haar measure of the pion fields, entering the path integral, generates 
an effective Lagrangian 8Ch which is generally singular in the continuum limit. There exists one 
parameterization of the pion fields only, for which the weight factor is equal to unity and SCh = 0, 
respectively. This unique parameterization ensures selfconsistency of the MF approximation. We 
use it to calculate thermal Green functions of the pion gas in the MF approximation as a power series 
over the temperature. The Borel transforms of thermal averages of a function J(x a X a ) 01 the pion 
fields x° with respect to the scalar pion density are found to be -j=Cf(4t). The perturbation series 
over the scalar pion density for basic characteristics of the pion matter such as the pion propagator, 
the pion optical potential, the scalar quark condensate < qq >, the in-medium pion decay constant 
F, and the equation of state of pion matter appear to be asymptotic ones. These series are summed 
up using the contour-improved Borel resummation method. The quark scalar condensate decreases 
smoothly until T max ~ 310 MeV. The temperature T max is the maximum temperature admissible 
for thermalized non-linear sigma model at zero pion chemical potentials. The estimate of T ma x is 
above the chemical freeze-out temperature T ~ 170 MeV at RHIC and above the phase transition 
to two-flavor quark matter T c ~ 175 MeV, predicted by lattice gauge theories. 

PACS numbers: ll.10.Wx, 13.25.Cq, 14.40.Aq 



I. INTRODUCTION 



Ultrarelativistic heavy-ion collisions allow to study phase diagram of QCD at high temperatures and small chemical 
potentials. Under such conditions, lattice gauge theories (LGTs) predict a crossover or a first-order phase transition 
into the deconfined phase, accompanied by the restoration of chiral symmetry. In two- and three-flavor LGTs, the 
critical temperature is estimated to be T c ~ 175 MeV and T c ~ 155 MeV with a 5% systematic error HHHQ. These 
values are close to the chemical freeze-out temperature determined from statistical models by fitting particle yields 
in heavy-ion collisions |4|. The chemical freeze-out temperatures T — 160 -j- 174 MeV and T — 160 -j- 166 MeV are 
extracted [| at top SPS energies (^/snn — 17.3 GeV) and top RHIC energies (y/spjN — 200 GeV), respectively. The 
boundary of the phase transition might apparently be already crossed. 

In ultrarelativistic heavy-ion collisions pions are the most abundant particles 0, 0, @ 1 so thermodynamic char- 
acteristics of a pion-dominated medium and in-medium pion properties are of current interest. At zero chemical 
potentials, light degrees if freedom, i.e. pions, give the dominant contribution to the pressure of hadron matter 0, 
while contributions of resonances like p- and w-mesons with masses M ^> T c are suppressed as exp(— ¥r) <C 1. 
According to the Gibbs' criterion, the balance of pressure determines the critical temperature of the phase transition 
at fixed chemical potentials. The QCD phase transition can be driven therefore by pions i.e. by the lightest QCD 
degrees of freedom. 

The energy density of a resonance gas is proportional to the small parameter exp(— ^) <C 1 which can, however, 
be compensated by an exponentially large amount of resonances, as has been conjectured by Hagedorn |9j, llOj and 
discussed recently e.g. in Ref. The contribution of resonances to the pressure is still suppressed in such models 
by a factor <C 1 where M ~ 1 GeV is a typical scale of resonance masses. 

Pions are Goldstone particles and play a special role in the restoration of chiral symmetry. LGTs give an evidence 
that deconfmcment and chiral phase transitions occur at the same critical temperature for zero chemical potentials 

m 

Chiral perturbation theory (ChPT) exploits the invariance of strong interactions under the SU(2)l <g) SU(2)r sym- 
metry [l 2L ll3( . It has been proven to be highly successful in phenomenological descriptions of low-energy dynamics of 
pseudoscalar mesons [l4l Il5| . The pion matter created in relativistic heavy-ion collisions is charge symmetric £jj , so 
the pion chemical potentials can be set equal to zero. The ChPT expansion in the vacuum runs over inverse powers of 
the pion decay constant F — 93 MeV and, in the medium, additionally over the pion density (or temperature). Prop- 
erties of the pion gas at finite temperatures, the associated chiral phase transition, and the in-medium modifications 
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of pseudoscalar mesons within ChPT have extensively been studied EH E3, EM H3, [53, H^, . 

The optical potential of an in-medium particle yields the self-energy operator to the first order in the density. It is 
determined by the forward two-body scattering amplitude on particles of the medium. ChPT is suited for calculation 
of the pion self-energy beyond the lowest order in density, since ChPT has all multi-pion scattering amplitudes fixed. 
The forward scattering amplitudes of a probing pion scattered off n thermalized pions describe 0(p n ) contributions 
to the in-medium pion self-energy operator. The leading ChPT amplitudes do not contain loops. This approximation 
is equivalent to the mean field (MF) approximation. 

In this paper, we perform to the leading ChPT order the generalized Borel resummation of the asymptotic density 
series for the pion propagator, the pion optical potential, the scalar quark condensate < qq >, the in- medium pion 
decay constant F, and the equation of state (EOS) of pion matter. 

The outline of the paper is as follows: In the next Sect., we find for the pion field a parameterization, which 
provides a constant Lagrange measure for path integrals. The integration over the pion field variables is extended 
from — oo to +00 to convert path integrals to the Gaussian form. Such a procedure does apparently not affect the 
perturbative part of the Green functions. In Sect. Ill, a method for calculation of thermal averages within ChPT 
in the MF approximation is described. The thermal averages appear as inverse Borel transforms with respect to the 
non-renormalized scalar pion density. The lowest order expansion coefficients are calculated for the pion effective 
mass and other quantities and compared to the earlier calculations. In Sect. IV, we show that the power series with 
respect to the density are asymptotic and hence divergent. They can, however, be summed up with the help of the 
contour-improved Borel resummation technique. In Sect. V, predictions of the non-linear sigma model are compared 
to LTGs. The results obtained by the resummation are discussed in Conclusion. 



II. PION FIELD PARAMETERIZATION 



The lowest order ChPT Lagrangian has the form 

F 2 F 2 M 2 

C = — Tr^Ud^} + -^Tr[U^ + U] (II.l) 

where M n is the pion mass. The kinetic term for the matrix U(x) G SU(2) is invariant under the chiral transformations 
U — > U' = RUL + where R, L £ SU(2). The second term in Ea. HII.ljl breaks chiral symmetry explicitly. In what 
follows, we set F = 1. 

The matrix U(x) can be parameterized in various ways. The on-shell vacuum amplitudes do not depend on the 
choice of variables for pions [27], . The method proposed by Gasser and Leutwyller E3 establishes a connection 
between QCD Green functions and amplitudes of the effective chiral Lagrangian. Using this method, the QCD on- 
and off-shell amplitudes can be calculated in a way independent on the parameterization. The detailed studies of Refs. 
[2^.l30| demonstrate that the in-medium effective meson masses are independent on the choice of meson field variables 
up to next-to-leading order in ChPT and to first order in density, in accordance with the equivalence theorem. The 
in-medium off-shell behavior and going beyond the linear-density approximation are discussed in Refs. 

In order to provide the S'-matrix invariant with respect to a symmetry group, both, the action functional and the 
Lagrange measure entering the path integral should be invariant. The chirally invariant measure dfi[U] — dfi[RUL + ] 
entering the path integral over the generalized coordinates coincides with the Haar measure of the SU(2) group. In 
the exponential parameterization, 

U(<t> a ) = e lr °* Q (II.2) 

and 

dfi[U] = sm 2 ((t))(f)~ 2 d 3 (f) (11.3) 

(see e.g. [33|). The leading order ChPT is equivalent to the 0(4) non- linear sigma model due to the isomorphism of 
algebras su(2)l su(2)r ~ so(4). The quantization of the 0(4) non-linear sigma model [34[ yields automatically the 
measure (|II.3|I . 

Any perturbation theory uses for dynamical fields an oscillator basis to convert path integrals into a Gaussian form. 
It is necessary to specify variables, in terms of which the Lagrange measure is "flat": 

dfJ,[U] = d 3 X - (II.4) 

The weight factor can always be exponentiated to generate an effective Lagrangian 5Ch, in which case \ a = 4> a 
provides the desired parameterization. The exponential parameterization (|II.2fl gives, in particular, SCh = 
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— log(sin 2 (0)0~ 2 ) where a is a lattice size. SCh diverges in the continuum limit. The non-linear sigma model 
is not a renormalizable theory, so divergences cannot be absorbed into a redefinition of F and M n . Using the MF 
approximation, it is usually possible to keep renormalizations finite. The exponentiation of a variable weight factor 
breaks, in general, sclfconsistency of the MF approximation. 

The divergences arising from SCh could, however, be compensated by divergences coming from the higher orders 
ChPT loops. From this point of view it looks naturally to attribute SCh to higher orders ChPT loop expansion 
starting from one loop. The MF approximation for the ChPT implies then that the tree level approximation neglects 
SCh from the start. It is hard to expect that such an approximation is relevant at the high temperature limit where 
the chiral invariance is supposed to be restored. 

The interaction terms in the effective Lagrangian which appear due to presence of the Haar measure have been 
discussed earlier in QCD |35L l3t| . The exponentiation of the weight factor gives consistent results due to renormal- 
izability of QCD. The divergences appearing in the continuum limit from SCh at a tree level are compensated by 
one- loop gluon self-interaction diagrams. 

The consistency of the MF approximation of the non-linear sigma model survives with one parameterization only, 
which uses the dilatated pion fields variables \ a = Q x/</ ) such that \ 2 x' — sin 2 (0) > where \ ~ (x Q X Q )^ 2 and 

X 3 = ~(<£-sm(<£)cos(</>)). (II.5) 

The vacuum value 4>" ac = corresponds to Xvac = 0, x is a monotonously increasing function of (f>. The value of 
Attx /3 has the meaning of a volume covered by a 3-dimensional surface of radius x. m a 4-dimensional space. 

The parameterization based on the dilatation of <j> a gives SCh = 0, does not require the higher orders ChPT loop 
expansion for the consistency, and allows to work in the continuum limit with finite quantities only. 

The SU(2) group has a finite group volume. The magnitude of the x° fields is restricted by Xmax = (37T) 1 / 3 . It is 
generally believed that using the perturbation theory, one can extend the integrals over x a from — oo to +oo. The 
modification of the result is connected to large fields fluctuations of a non-perturbative nature, which do, apparently, 
not affect perturbation series. This conjecture is essential for the standard loop expansion in ChPT. A similar 
extension of the integration region is used in QCD [35l l36| . 

The Weinberg's parameterization |37| does not restrict the magnitude of the pion fields. Such a parameterization 
looks especially attractive as it simplifies conversion of the path integrals into the Gaussian form. Because of the 
variable Lagrange measure, it can be effective starting from one loop. 

The propaga tors (heat kernels) of free particles moving on compact group manifolds have been analyzed in Refs. |38l 
IME31II1E1 It was shown that the semiclassical approximation for the propagators is exact. An explicit form 
of the propagator for the group SU(2) is given by Schulman |^£j and Duru 0. The problem of the variable 
Lagrange measure of the path integrals appears already at the quantum mechanical level. The path integral as it 
was recognized by Marinov and Terent'ev |3^| is ill posed at the tree level. As demonstrated by Baaquie jl^, the 
divergent contributions coming from the variable Lagrange measure are cancelled by loops generated by the effective 
Lagrangian. A system of coupled oscillators on a compact group manifold represents unsolved problem. 

The requirement of chiral invariance of the Lagrange measure dfi[U] being combined with the requirements that 
(i) the perturbation expansion is based on the oscillator basis and (ii) the MF approximation does not involve loops 
restricts the parameterizations of the pion fields to only one admissible parameterization. 

III. THERMAL AVERAGES IN THE MF APPROXIMATION 

The quadratic part of the Lagrangian with respect to the fields x" can be extracted from Eq. ijll.lfl and the rest SC 
treated as a perturbation: 

£ = 2 W - -fx 2 + SC (III.l) 

where M w is an effective pion mass to be determined self-consistently and 

1 1 

S£ = £ W Jo + ^X^xVxV 2 Ji + Ml(J 2 - 1) + -fj 3 , (III.2) 
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with 



Jo = ^-l, (IH.3) 

X 



^--4t-t-^' ( IIL4 ) 



X 4 _ sin 2 ( 0) 
sin 4 (0) X 2 

J 2 = cos(0), (HI.5) 

Jz = x 2 - (in-6) 

The physical observables are expressed in terms of the thermal Green functions. 



A. Thermal Green functions 

The scalar pion density p in a thermal bath with temperature T is defined by 



<>->=•• W .-='>-' A / 7^a2n (k) (III.7) 



where Z x is a renormalization constant of the pion x-field and 



The average <|III.7ll involves the normally ordered % a operators. The energy w(k) = y M. 2 . + k 2 which enters the 

Bose-Einstein distribution contains the in-medium pion mass M w . A similar situation occurs in the Bogoliubov model 
of a weakly interacting non-ideal Bose gas. Also in Fermi liquid theory the momentum space distribution is determined 
by the in- medium dispersion law of quasi-particles (see e.g. 
In terms of thermal two-body Green function, 

g^(r, x) = - < T X a (r, x) X ' 3 (0, 0) >, (III.9) 

the momentum space distribution noQs.) is determined by equation 

r%(2n (k) + ^) = -T £ lim o £ Q^{u> s , k) e "-. e (111.10) 



s — — oo 



where w, = 2ttsT and 



0^(r,x)=T^ / __0^(o; s ,k)e-^ + * x . (III.ll) 
The density p can be found from 

5 a ^(2p + ,w) - - lim g^(e, 0), (111.12) 

e^ — 

where /9,j ac is the vacuum density of the zero-point field fluctuations 

dk 



Pvac — Z 



2w(k)(27r) 3 ' 



The vacuum density can be treated self-consistently beyond the lowest order ChPT only. In the MF approximation 
it is assumed that zero-point field fluctuations are absorbed to a redefinition of physical parameters entering the 
Lagrangian. In our case, the pion mass and the pion decay constant receive these contributions. 

We thus neglect the vacuum fluctuations. In what follows, p vac is set equal to zero wherever it appears. This is 
equivalent to using thermal part of the Green function within the loops: 

lim g a/3 (e,0)^ lim g a(3 (e,0)- lim g a(3 (e, 0)| T=0 . (111.13) 

e^ — e— > — e— ^ — 

Ea. i|III.12|l is illustrated graphically on Fig. ^ Since we neglect the zero-point field fluctuations, we may not worry 
on the ordering of operators, entering thermal averages such as 1III.7L (IIII.14II and others, anymore. 
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FIG. 1: Diagram representation of the thermal average < X a X t> > defined by Eqs. 1111,711 and 1111. 1211 . The thick dashed line 
represents the dressed pion propagator. 



B. In-medium pion mass and other observables 

In the MF approximation, the self-energy operator can be calculated from 

S a ^{x -y) = -< - f - / SC(z)d 4 z > . (IIL14) 

Equation < x a ®nX^ > = is a consequence of the symmetry of the pion matter with respect to isospin rotations. In 
the momentum representation, using equation 



one gets 



where k 2 = —to 1 — k 2 and 



< d^ X a d^x f3 >= S a ^M 2 2p, (111.15) 
£(fc 2 ) = fc 2 S! + A/ 2 £ 2 + Af 2 £ 3 (HI.16) 



51 = - < Jo + \ji >, (HI.17) 

5 2 = ^-<(l-^)J 2 >, (HI.18) 

2p bp 

£ 3 = \ < (1 - ^)G7o + \jx) > -1- (HI.19) 



The self-energy operator is presented graphically on Fig. [3 



>' + W +■ 



FIG. 2: Diagram representation of the self-energy operator in the MF approximation. 

Eg . (jlll. 15|) can be derived as follows: The thermal average (|III.15ll can be rewritten in terms of the two-body Green 
function 

< 9,/^/ >= lim □e Q/3 (r,0) (111.20) 

where □ = —^3 — A. The two-body Green function satisfies equation 

(□ + A/ 2 + £(-D))£ Q ^,x) = -rt(r)«5 3 (x).. (111.21) 

The perturbation 8L(z) is quadratic with respect to the derivatives, so X(fc 2 ) in the MF approximation is a first order 
polynomial. The expansion l|III.16|) makes this feature explicit. Ea. (|III.21(l can be rewritten in the form 

(□ + Ml)G af \T, x) = -Z x 5 al3 6(T)6 3 (x) (111.22) 

where Z x is the renormalization constant introduced earlier. Using lim T ^ +0 <5(r) = and Eas. ljIII.12p . (fill. 201 . 
(|III.22fl . and i|A.4|) . we arrive at Ea. ljIII.15fl with 2p replaced by 2p + p vac . The value of p vac must be neglected, 
however. 



The effective pion mass can be determined from equation S(M^) = 0: 



Ml = -M* ' , (111.23) 

^1 + ^3 



while the renormalization constant is given by 



x dk 

The Green function in the MF approximation has the form 



Z;' : I-^1P = 1-S, (111.24) 



G a0 (Lu s , k) = -<5 Q/3 ^ ^. (111.25) 



It is presented graphically on Fig. [3J 



FIG. 3: Diagram representation of the dressed pion propagator (thick dashed lines) in terms of the scalar density and the bare 
pion propagators (thin dashed lines) . The scalar density appears from loops formed by the dressed pion propagator according 
to Ea. ftH.12l . 

In what follows, we need also functions 

J a = — 111.26 

X 

J 5 = -J2(-J4 + (20p-2 X 2 )j; + 8p X 2 J") 

+J 4 ((12p - 2 X 2 ) j£ + 8p X 2 J 2 ") (111-27) 

which appear in calculations of density dependent renormalization constant Z„ and the pion decay constant F ac- 
cording to Eqs. (|III.29|> and (|III.31(1 . Here, J[ and 3" are first and second derivatives with respect to <j). 

The thermal Green function of the pion is defined by & aP {x - y) = - < Tn a (x)^ (y) >. The pion fields ir a (x) = 
^Tr[r a U(x)] do not depend on derivatives of x° ■ I n such a case, no additional k 2 dependence appears as compared 
to the x propagator. Given the renormalization constant Z x , the renormalization constant Z n of the pion propagator 
can be found from 

5 a ^Z 1 J 2 Z- 1 ' 2 =< ^ > . (111.28) 



The pion propagator A a "(x — y) is shown graphically on Fig. 




FIG. 4: Diagram representation of the pion propagator A a ^(x — y). The dressed pion propagators are shown as thick dashed 
lines and the bare pion propagators are shown as thin dashed lines. The second line contains insertions into the ends of the 
dressed pion propagators only. 

The quark scalar condensate < qq > is proportional to the expectation value of the scalar field a = hTr[U(x)]. The 
calculation of the averages < d~K a jdx > and < a > leads after factorization of isotopic indices to the calculation of 
averages defined in terms of functions 1)111.5(1 and ((111.26)1 : 

Zl^Z- 1 ' 2 = < ^-Ji >, (111.29) 
A bp 

<qq> I <qq > vac = < Jl > ■ (111.30) 
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<x > 



\a a' 



, a |5\ 



FIG. 5: Diagram representation of the thermal average < \ A >=< X a X°'X^X^ >~ (3x3 + 3 + 3) x 4p 2 . Greek indices e.g. a 
and (3 within one loop designate S a/3 . The result 60p 2 is in agreement with Eg. 1111.321 for n — 2. 



In the pion gas, the spectral density of the axial two-point function < T A ( ^(x)A l l(y) > is characterized by two 
different pion decay constants |23| |. In the limit T = 0, these constants coincide with F. Let us find an in-medium 
pion decay constant, F, appearing in the spectral density of the two-point function < Td fl A"(x)d l/ A^,(y) >.. This 

definition is equivalent to equation d^A^ = FM^Z^ 1 ^ 2 n a where the right side gives the physical pion Z w 1 ' 2 n a .. The 
axial current has the form A^ = —ad^ir" + (d tJl a)TT a . The calculation of the average < TA^(x)A^(y) > leads after 
factorization of isotopic indices to the calculation of an average defined in terms of the function i|III.27|) . In this way 
we obtain 

FZ- 1 ' 2 =< ^-J 5 > . (111.31) 



C. Thermal averages as inverse Borel transforms 

Suppose we want to find the thermal average of a function J(x 2 ) which can be expanded in a power series of \ 2 
The average values of powers of the operator x 2 = X 1 X 1 can be calculated as follows: 

' si!s 2 !s 3 ! 

sis 2 s 3 

= E -r ! V-l(2 S i-l)H(2 S2 -l)!!(2 S 3-l)!!(2 /9 ) n 

S1S2S3 



1 1 r°° r 00 c 00 111 

E -rVr4^( 4 P)" / / / xr^x 8 r^x 8 r^e-^^dx 1 dx 2 dx, 



™ sl!s2!s3l7r3/2 ' ■'" 



' JO JO JO V^2^3 



n 



4 

A 
2 

A Jo 



(4px 2 )"a; 2 e- a; 

(ipt^t^e^dt. (111.32) 



The summations run for si + S2 + S3 = ^, where s, = 0, 1, .... In the first line, we use the multibinomial formula for 
((x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 )™- The factorization of the average in the second line is a consequence of the MF approximation 
according to which the averages are calculated over the free gas of quasiparticles, i.e., non-interacting effective pions 
in three isotopic states. In order to calculate, e.g., the average value < (x x ) 2si >> one has to consider (2si)! 
permutations of the operators x 1 due to their possible pairings. Permutations inside of each pair are counted twice, 
whereas permutations of the pairs are counted S\\ times. This gives the factor (2si)!/(2 Sl Si!) = (2s\ — 1)!! in the 
third line of 1III.32|) . To arrive at the final result one has to express (2sj — 1)!! in terms of Euler's gamma functions 
and use their integral representations to convert tree-dimensional integral into the one-dimensional integral. 

The MF approximation we used consists in the neglection of dynamical loops, i.e., loops composed from more than 
one dressed pion propagator. The loops composed from one dressed pion propagator give the density (|III.7|1 after 
neglecting p vac . The calculation of the n — 1 average value is graphically presented on Fig.. ^ (the indices a and 
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/3 have to be contracted with 8 a @). The calculation of the n = 2 is illustrated on Fig. |S| An alternative proof of 
Eci. (|111.32l) . based on the thermal Green functions method, is given in Appendix A. 
The average value of J(x 2 ) can be written in the form 

< J(X 2 ) >= -t= / J(4pt)tV 7 e- t dt. (111.33) 
V 71 " Jo 

The thermal fluctuations are suppressed exponentially according to the value of x 2 ( see a ls° Ea. ljB.ljl V For small 
fluctuations, x coincides with (f>. 

The thermal average < J(x 2 ) > to all orders in the pion density and to leading order in ChPT is given by the 
inverse Borel transform of the function -j=J{b£). 

Using the explicit form of J7i and re-expanding results in terms of the physical density po = Z~ 1 p 1 we obtain for 
Po « 1 

„V2/»r 2 , 39 , 2601 o 614163 , , TXT „^ 

Ml/Ml = 1 + po - -pi - —pi - -j^-pt + (IIL34) 

, 12 2 368 , 80832 , 

Z- 1 = 1 + Y Po + ^5 Po + -^5~ Po + -' (IIL35) 
„ 1 „ 46 2 1108 o 204202 4 

= 1 + 2 P0 + + __ p 3 + __ p 4 + ^ (nI36) 

/ - 1 o 3 2 39 o 7803 4 , nim , 

< qq > I < qq > vac = 1 - ip Q - -p 2 + —p* + ~^Pl + (HL37) 

F/F = 1 - 2 P0 + \pl + ^§ P l + + ... • (in-38) 

For massless pions po = j^T 2 . In such a case, the power series expansions over density and temperature coincide. 

The result for M 2 to order O{po) is in agreement with 0, H^l, |2£j , the result for Z^ 1 to order O(po) is in agreement 
with 25], the result for < qq > to order 0{p\\) is in agreement with 0,0,0,0,113, the result for F to order O(po) 
is in agreement with (Hill 1^. 

Chiral invariance of the non- linear sigma model in the MF approximation is commented in Appendix B. 



IV. SUMMATION OF THE DENSITY SERIES 

Most perturbation series in quantum field theory are believed to be asymptotic and hence divergent 0, 0] . The 
integral representation l|III.33|l permits by changing the variable t — ► t' = t/p an analytical continuation into the 
complex half-plane e[p] > 0, whereas at e[p] < the integral diverges. The half-plane e[p] < contains in general 
singularities the character of which depends on properties of the function J(x 2 )- The convergence radius of the Taylor 
expansion is determined by the nearest singularity. One cannot exclude that the point p = around which we make 
the expansion is a singular point. 

Inspecting Eas. l|III.3|l - (|III.6|l and l|III.26|l - (|III.27Jl . we observe that the Borel transforms involving J\ and J§ are 
singular at 4> = kit where k = 1,2, .... The integral along the real half-axis of the Borel variable t is therefore not 
feasible. Series involving such functions are asymptotic, have zero convergence radii, and are furthermore not Borel 
summable [H0|47|. 

In order to clarify the character of the power series over the physical density po, we calculate the higher order 
expansion coefficients. We write Ojj °o) = Sfe°=o c nPo f° r observables (|III.34|I -I |III.38|) and plot on Fig. Elthe ratios 
c-n+i/cn versus n up to n — 50. 60] The apparent asymptotic regime starts at n ~ 10. The ratios c„+i/c„ increase 
linearly with n, indicating clearly that the power series are asymptotic. The slopes of these ratios are approximately 
equal. One can expect that the nearest singularity <j) = ir that appears in the series expansion over p appears again in 
the series expansion over po by virtue of p = Z x po. In such a case, we would expect a singularity in the Borel plane 
at U = (^f ) 2/3 as well. Fr om the other side, for c„ ~ n\a n the Borel transform is singular at t — 1/a. The slope 
equals then a — ) 2 / 3 ~ 1.4, in good agreement with results presented on Fig. |SJ The power series (|III.34|) - HIII.38() 
with respect to the physical density p$ have zero convergence radii and are also not Borel summable. 

Generalizations of the Borel summation method have been proposed and effectively used in mathematics and atomic 
physics 0,E3- In our case, fortunately, the residues &t <j) = kn of the Borel transforms vanish. A contour- improved 
Borel resummation technique can therefore be applied which consists in a shift of the integration contour into the 
complex t-plane. The result is stable against small variations of the contour around the positive real axis. The thermal 
averages, in particular, do not acquire imaginary parts and thus the expectation value, e.g., of the hermitian operator 
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FIG. 6: The ratios c n +i/c n of the expansion coefficients of power series X]^Lo c ™Po over the pion density po for (a) the in- 
medium pion mass M%, (b) the propagator renormalization constants Z~ x of the x-field and Z~ x of the pion field, (c) the scalar 
quark condensate < qq > , and (d) the pion decay constant F. The linear growth of the ratios beyond n ~ 10 indicates that 
all power series are asymptotic. The same sign of c n for n > 10 indicates furthermore that the series are not Borel summable. 
The approximately equal slopes are related to a singularity of the Borel transforms on the real axis at 4> = T (see text). 



qq remains real. The uniqueness of the result justifies the integrations by parts of the functions performed to 
derive Eos. (1111. 1811 . 1111.1911 , 1111.2911 . and 11111.3 H . 

The power series <|III.34|) - (|III.38|) become summable using the contour- improved Borel resummation method. 

The EOS of pion matter can be determined by averaging the energy-momentum tensor constructed from the 
Lagrangian (|III.lJl . The renormalized energy density e = Z~ x Tqq and the pressure p = Z~ x J2i=i 3 Tu/3 are given by 

£ = 6 / (2?i^ 2(k) " o(k) " Z * * ( AI * 3p0 + M " {<T ~ 1} ) ' (IV ' 1} 
P = 2 J J^^M^+Z-^M^po + M^a-l)) (IV.2) 

where a =< qq > / < qq > vac - The results are plotted in Fig. 

There exists phenomenologically a broad interval of < T < 2.48F ~ 230 MeV with Z x < 1, covering the range 
of temperatures specific for heavy-ion colliders. Above T ~ 230 MeV the renormalization constant Z~ x < 1. In 
the vacuum, one has Z^ 1 > 1 as a consequence of positive definiteness of the spectral density of two-point Green 
functions. In the medium, the spectral density of temperature Green functions of bosons is not positive definite (see 
e.g. HI, Chap. 17-3), so Z x 1 < 1 does apparently not contradict unitarity. Although the product of uj and the pion 
spectral density cannot be negative, sum rule for this product is unknown. The region of self-consistency of the MF 
approximation within perturbation theory over the density extends up to T max — 3.36-F ~ 310 MeV. 

An interesting finite-temperature extension of the Weinberg sum rule for the product of lo and the spectral den sity 
of vector mesons is discussed in In the MF approximation resonances do not exist, so the constraints 

apply starting from one loop ChPT. 
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FIG. 7: The non-renormalized scalar pion density p, the renormalized scalar pion density po, the ratio M-njM^ between the 
in-medium and the vacuum pion mass, the inverse renormalization constants Z^ 1 and Z^ 1 of the pion x-held and the pion 
field ir a (x) = ^Tr[T a U(x)], respectively, the scalar quark condensate < qq >, and the in-medium pion decay constant F versus 
temperature T are shown. The energy density e (dashed curve) and the pressure p (dot-dashed curve) of the pion gas are also 
plotted. Solutions do not exist above T max = 3.36F ~ 310 MeV (the right vertical line). The critical temperature T c ~ 175 
MeV of phase transition into the quark matter is shown by the left vertical line. All quantities are given in units of the pion 
decay constant F — 93 MeV. 



For T < T max = 3.36.F ~ 310 MeV the system has regular behavior. The scalar quark condensate decreases 
smoothly to zero. The effective pion mass increases with temperature in agreement with the partial restoration of 
chiral symmetry and approaches M^^nax = 1.61M W . The pion optical potential V pt = {M% — M%) / '(2M ff ) increases 
also. 

The pion optical potential V op t is calculated to all orders in density. It corresponds to the summation of the forward 
scattering amplitudes of a probing pion scattered on n surrounding pions, with n running from one to infinity. 

For T ~ T, nax , Z~ l — > oo. The pions disappear from the spectral density of the propagator < T n a {x)^ (y) >. 
As elementary excitations, the pions do not disappear, however, from energy spectrum Ea. (|IV.l|l and contribute to 
the pressure Ea. (|IV.2ll . Moreover, the pions contribute to the spectral density, e.g., of the two-point Green function 
<Tir a (x)xHy) >■ 

The estimate of T max is significantly above the chemical freeze-out temperature T ~ 170 MeV at RHIC. 

It is worthwhile to notice that for massless pions, the scalar quark condensate vanishes at T max (M„ = 0) = 2.46^ ~ 
230 MeV. The behavior of observables 1|III.34|) -( |III.38(1 does not change qualitatively. EOS becomes identical with 
EOS of the ideal pion gas, as can be seen from Eas. (|IV.l|l and (|IV.2|I . 




FIG. 8: The normalized scalar quark density a =< qq > / < qq > vac versus the pion scalar density po- The power series 
expansion (II1I.37I is truncated to orders O(po) for n — 1,2,3 and 4 (long-dashed, dashed, dotted, and dot-dashed curves, 
respectively) and results are compared to the exact numerical summation of the series HIII.37(I using the contour-improved 
Borel resummation technique (solid curve). The value of po, c is the critical scalar pion density for phase transition into the 
quark matter. The value of po,max is the maximum scalar pion density admissible in thermalized non-linear sigma model. 



Let us check accuracy of presentation of the results by truncated series. In Fig. we show the ratio a =< qq > / < 
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qq >vac calculated using the series expansion (|III.37|I truncated to orders O(pq) for n = 1,2,3 and 4. It is compared 
to the exact (numerical) summation of the perturbation series. We see that in the hadron phase T < T c , the first 
order approximation is very precise. Noticeable, irregular deviations appear beyond T c only. We expect, therefore, a 
20% decrease of the scalar quark condensate at T — T c . 

The first-order approximation is sufficiently precise for M% also: At T = T c , the first order of l|III.34|) gives for 
the effective pion mass M% ~ 1.07M^, while the exact summation gives 1.03M^. The other observables display the 
similar features. 

A prescription for approximate summation of asymptotic series, which is effective if initial sequential terms of 
asymptotic series decrease first before increasing, can be found in Ref. [50|. One should attempt to truncate asymptotic 
series where two sequential terms are of the same order. The magnitude of the first neglected term gives the accuracy 
of the summation. Comparing n = 1 and n = 2 terms of <|III.37|) . one may conclude that accuracy of the series 
truncated at n = 1 is 100% for po — 2F 2 . One can expect that for po c — 0.07F 2 <C 2F 2 , the accuracy is very good. 
As we observed, this is the case. The same arguments require po be less than 0.25F 2 for M 2 . Such a requirement is 
satisfied also, however, with a lower precision. 



V. COMPARISON WITH LATTICE GAUGE THEORIES 



LGTs allow to calculate QCD observables from first principles. The smallness of physical quark masses and finite 
lattice spacing restrict the power of LGTs to temperatures above ~ 100 MeV. The non-linear sigma model, on the 
other hand, is an effective theory of QCD at temperatures small compared to the pion mass. There are no intrinsic 
restrictions to extrapolate results of the non-linear sigma model up to T max ~ 310 MeV. However, above pion mass 
the non-linear sigma model represents in the strict sense a model rather than an effective theory of QCD. 

The phase transition to the quark matter appears in LGTs at T c ~ 175 MeV PIQ- The domain of validity of the 
non-linear sigma model is therefore limited to temperatures below T c . The interval of temperatures from ^100 MeV 
up to T c is suitable for a meaningful comparison of predictions from the non- linear sigma model and LTGs. 




FIG. 9: The scalar quark condensate a =< qq > / < qq > vac versus the temperature in units of F = 93 MeV obtained within 
the non-linear sigma model using the MF approximation (solid line below T c and dashed line above T c ) and from lattice gauge 
theory for two flavors [Tl|. 



Fig. El shows the temperature dependence of the quark scalar condensate. The overall mass scale in LGTs is fixed 
assuming the string tension is flavor and quark mass independent. The agreement of the non-linear sigma model 
and LTG [TlT | is not unreasonable with regard to first two-three points with lowest temperatures. In the narrow 
deconfinement region the models strongly diverge, basically because the non-linear sigma model does not expose 
quark-gluon degrees of freedom. 

The energy density as a function of temperature is shown in Fig. ^|at zero chemical potentials. The solid line and 
the dashed line are predictions of the non-linear sigma model. In addition, the lattice simulations from Refs. [n| 
are shown. 

At a first-order phase transition, the energy density experiences a jump while the pressure remains a smooth 
function. It is not quite clear from the LGT data presented on Fig. EH whether we observe a crossover or a first-order 
phase transition. Usually, in the two-flavor case lattice simulations give evidence for a first-order phase transition, 
while pure gauge theory predicts a second-order transition and 2+1 flavor LGTs with physical quark masses predict 
a smooth crossover 13. 
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FIG. 10: The normalized energy density e/T 4 versus the temperature T obtained within the MF approximation of the non- 
linear sigma model (solid line below T c and dashed line above T c ) and from lattice gauge theory (filled triangles) with two 
flavors [H ITl|. The temperature is given in units of F = 93 MeV. 



If the phase transition is of first order , the energy density in LTG appears to be at T c several times greater than in 
the non-linear sigma model. Such a divergence can be attributed to resonances like p- and w-mesons which contribute 
to the energy density, but do not exist in the non-linear sigma model in MF approximation. An ideal resonance gas 
model which is in agreement with lattice simulations has e.g. been developed by Karsch et al. 

In the case of a first-order phase transition, the jump in the energy density can be as large as As c /T 4 ~ 8, in which 
case pions alone saturate the energy density. If the jump is smaller, Hagedorn's conjecture on the exponential growth 
of the number of resonances with mass might be appropriate. The value of the jump is crucial in order to understand 
the role of higher resonances. 

Modelling the first-order phase transition within the framework of the MIT bag model jHJ allows to determine 
the critical temperature using the pressure balance pq = pu + B where B — 57 MeV/Fm 3 is the vacuum pressure 
in 2 + 1 flavor QCD |52^. pq and pn are pressures in quark and hadron phases, respectively, and determine jump in 
the energy density Ae c = eq + B — er where eq and Eh are energy densities of quark matter and hadron matter, 
respectively. For ideal gases of relativistic particles, Ae e = 3po + B — 3pn = 42? , and so Ae c /T 4 ~ 3. The vacuum 
pressure B is density and temperature dependent IHbI IfiH l57l IHflj . which brings to Ae c more uncertainties. 

The energy density predicted by LTGs is obviously underestimated around T ~ 100 MeV due to unphysical quark 
masses mq/T = 0.4 and restrictions from finite lattice size. The energy density and the pressure depend strongly on 
the quark masses and the pion mass. The horizontal part of EOS at high T (dashed line and partially solid line) is 
close to the energy density of massless pions e/T 4 — 7r 2 /10. 









NLSM 






LGT 






hadron matter 


/ quark matter 











T/F 

FIG. 11: The normalized pressure p/T 4 versus the temperature T in units of F = 93 MeV in the non-linear sigma model 
(NLSM) using the MF approximation (solid line below T c and dashed line above T c ) and in lattice gauge theory (LGT) with 
two flavors [ltlTl|. 



The pressure is plotted on Fig. ^2 The lattice predictions 0, [n| and the non-linear sigma model predictions are 
in agreement at the phase transition point, while at smaller temperatures the pressure in LGTs is significantly lower. 
The pressure of the pion gas approaches the pressure of ideal gas of massless pions p/T 4 = 7r 2 /30 at T — > T max . 
However, it does not fully reach the ultrarelativistic limit, so the effects from the finite pion mass remain. 
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According to the Gibbs' criterion (see e.g. J5JJ) the phase of matter with the highest presssure at equal temperatures 
and chemical potentials is preferable. At T ~ T c , heavy resonances are still slowly moving particles which contribute 
to the energy density Se mainly through their rest mass M and contribute to the pressure as Sp ~ j^Ss <C Se. For 
zero chemical potentials, the pressure in hadron phase is dominated by light pions [7j. 

The non-linear sigma model predicts the critical pressure p c reasonably well. The critical temperature can be 
obtained phenomenologically from Fig. 1111 applying the Gibbs' criterion for the hadronic EOS derived using the 
non-linear sigma model and the quark matter EOS derived from LGTs. The intersection of the two pressure curves 
yields then the temperature of the phase transition. The value obtained in such a way is in the remarkable agreement 
with lattice estimates based on the position of the steep rise of the energy density seen on Fig. ^| 



VI. CONCLUSION 



It is generally believed that ChPT represents an adequate tool for studying the pion matter at low temperatures. 
The Haar measure appearing in the path integral is important to keep the crural symmetry at high temperatures 
unbroken within the MF approximation. In an arbitrary parameterization of the pion fields one faces a dilemma: 
Accounting for the effective potential SCh arising due to the exponentiating the weigh factor of the path integral 
measure brings divergences which can be compensated by going beyond the MF approximation only, i.e., by including 
pion loops appearing in the higher order ChPT expansion. If SCh is neglected, the MF approximation does not 
restore the chiral invariance with increasing the temperature. There is only one parameterization with SCh = 0, 
which makes the MF approximation selfconsistent: The renormalizations are finite and the chiral symmetry is restored 
at high temperatures. This parameterization was used to study the in- medium modifications of pions and collective 
characteristics of the thermalized pion matter in the MF approximation. 

We made essentially two approximations: 

(a) The results are obtained by extending the integration region over the pion fields \ a from — oo to +00. 

(b) The MF approximation consists in the neglection of loops composed from more than one dressed pion propagator. 
Loops formed by one dressed pion propagator give the pion density (|III.7I) . i.e., the thermal parts of loops are accounted 
for, whereas the vacuum contributions of the zero-point field fluctuations are systematically neglected. An accounting 
for the vacuum parts of the pion loops in a thermal bath is possible beyond the lowest ChPT order. 

We constructed Borel transforms of most important thermal averages. The corresponding power series with respect 
to the density are asymptotic, have zero convergence radius, but are summable using the contour-improved Borel 
resummation method.. The pion gas does not exist above T max ~ 310 MeV, whereas at T < T max thermodynamic 
observables are smooth functions. 

The deconfinement temperature T c = 1.86F and the corresponding scalar density po,c = 0.07F 2 appear to be small 
enough for approximate however quite accurate calculation of observables in hadron phase using truncated asymptotic 
series. 

The region of validity of the method can be evaluated by requiring < 'X 2 >< Xmax- The corresponding restriction 
T < ^8(37r) 2 / 3 F w 560 MeV does not appear to be stringent. 

The method of summation of the density series proposed in this work can be extended to higher orders of ChPT 
loop expansion. 
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APPENDIX A: MF APPROXIMATION FOR < (x'x 1 T > 

The order of the operators \ a — X Q (0, 0) entering the thermal average < (x 7 X 7 )™ > is not specified so far. As 
we shall see, it is not important. Let us rewrite < (x 7 X 7 )™ > with the help of the ordering operator T . To make 
its action well defined, we set first arguments of the operators x"(0, 0) equal to small parameters e a , i.e., we make 
replacements x Q (0,0) — > x a ( e a,0). If the limit e a — ► exists and does not depend on the way the parameters e a 
approach zero, the thermal average < (x 7 X 7 )™ > can De expressed in terms of the Green function: 

< (X 7 X 7 )" >= I™ < ^X Q1 X Q1 X Q2 X" 2 -X Q "X Q " > • (A.l) 
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In order to evaluate the Green function, we pass to the interaction representation, 

<T X ai X ai X a2 X a2 -X an X an >=<Txo 1 Xo 1 Xo 2 Xo 2 -Xo"Xo"5(i 0) >< S(± 0) >~ x (A.2) 

where Xo( r i x ) = S(t,0)x° : (0,x.)S~ 1 (t,0) and S(t2,t\) is the thermal 5-matrix (see e.g. 0]). Applying the Wick's 
theorem to the 2n-body Green function (|A.2[i and neglecting all the diagrams except for those corresponding to the 
non-interacting gas of the dressed pions, we obtain 

S(ei,e 2 ,..,e 2n )= £ < Tx* X 71 >< T^ 2 X 72 > . . . < TV 3 " x 7 " > (A.3) 

(/3i/3 2 ...,3„),(7 l72 ... 7 „) 

where (/3i/?2- ■■/?«) and (7i72-.-7n) are permutations of (aiot2---ot n ). The summation runs over those permutations 
which occur according to the diagram decomposition of the 2n-body Green function of a non- interacting gas. The 
pion propagators are dressed. Entering Ea. i|A.3(l are therefore the Green functions < Tx°X 13 > instead of < T\q Xo >• 




FIG. 12: In the limit e a — > 0, the disconnected diagrams of the four-body Green function produce no loops with more than 
one dressed pion propagator, being thus significant in the MF approximation. 



Show on Fig. ^] are diagrams of the four-body Green function corresponding to the non-interacting gas of the 
effective pions. Such diagrams contribute to thermal averages. The diagram shown on Fig. 1131 which accounts for 
the pion rescattering should be discarded in the MF approximation. This feature holds for 2n-body Green functions. 
The MF approximation does not involve terms except for those entering Eq. ljA.31) . 




FIG. 13: In the limit e a — > 0, connected diagrams of the four-body Green function produce loops involving more than one 
pion propagator and therefore are not significant in the MF approximation. The only diagrams essential for thermal average 
< (x 7 X 7 )" > are those which correspond to the non-interacting gas of the dressed pions. 



The isotopic symmetry of the pion matter amounts to the two-body Green function Q a ^{r, k) which is symmetric 
under the permutation a <-> j3 at r = 0. Moreover, 

lim g al3 (e.k)= lim G al3 (e,k). (A.4) 

e->-0 

Eas. ljA.3fl and l|A.4|) imply that after neglection by all terms proportional to p vac the limit e a — > in Ea. (|A.l|l exists 
indeed and, at least within the MF approximation, does not depend on the way the parameters e a approach zero. 
Loops formed by one pion propagator give the density l|III.7|) . while loops involving more than one dressed pion 
propagator are neglected. 
We obtain therefore 

< (X 7 X 7 )™ >= (2p) n Yl c^ 171 ^ 272 ...^" 7 ". (A.5) 

(0 i/32.../3„),(7i72..-7n) 

The combinatorial structure of the Wick's decomposition is identical with that of the corresponding Gaussian integral, 
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so one can write 



<5 /3l71 (5 /32 72 ...(S /3 ' iT 

(/3l/?2— /3n),(7l72— 7 n) 

1 /• f f „ „ 1 



a;" 1 a;" 1 a;" 2 a;" 2 ... a;"™ a;"™ exp(--((x i ) 2 + (a; 2 ) 2 + (x'YVdx 1 dx J dx* 

(2TT)' i '' ! J J J 2 
4_ r+oo -i 

r 2n+2 exp( — r 2 )dr. (A.6) 



~ (2n)V 2 J ^ 2 

Replacing r — > V^t, we arrive at Ea. HIII.32ll 



APPENDIX B: CHIRAL INVARIANCE IN THE MF APPROXIMATION 



As discussed in Sect. II, the MF approximation imposes restrictions on the pion field parameterizations. If the tree 
level SCh is included into a MF calculation, one gets divergences. From other hand, the neglection by 8Ch violates 
the chiral symmetry. Let us demonstrate it explicitly: 

Equation (|III.33|I can be rewritten in more physical terms 

< >= / J(x 2 )e-- 2 ^ X . (B.l) 

The scalar product 

cos(6) = iTr[[/(0)[/t(0')] (B.2) 

is obviously chirally invariant. The value of = Q(<f>,<j)') determines the angular distance between two sets of the 
pion fields. The parameter space is thus a metric space with an infinitesimal distance 

dQ 2 = d(f> 2 + sin 2 0)(d(9 2 + sin 2 (9)dcp 2 ) (B.3) 

where 8 and ip are polar and azimuthal angles of the vector <f) a . 

Equation (|B.1|) has a transparent physical meaning. Quantum fluctuations of the fields x° induce fluctuations of a 
function J(x 2 )- These fluctuations are suppressed by the exponential factor entering l|B.l|) . 

According to Ea. i|II.5|) x depends on <f> which has the meaning of a metric distance 

^ = e(<r,cj (b.4) 

between the vector <p a and the vector </>" ac = which specifies the vacuum state on the 4-dimensional sphere, making 
thereby the chiral symmetry spontaneously broken. The only essential place in Eq. (|B.1|I where the dependence on the 
4>vac = shows up is the exponential factor. In the limit T — > oo the density increases, p — > oo, and so the dependence 
on 4>% ac = drops out. The exponential factor becomes a constant. It means that quantum fluctuations at all points 
of the 4-dimensional sphere contribute equally to thermal averages. This is in agreement with our expectations that 
the chiral symmetry restores with increasing the temperature. 

A use of the exponential parameterization with the Haar measure neglected would result in 

< ^) > J J( X 2 )e-^d% (B.5) 

The Euclidean volume d 3 <fi keeps explicitly a reference to the vacuum vector (j>y ac = 0. It means that the high 
temperature limit does not restore the chiral invariance in the MF approximation if the Haar measure is neglected. 
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